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ABSTRACT 

Aims. Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and 
when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (E(SFR)/M) 
in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift z « 1.6. 

Methods. We use deep far-infrared PACS data at 100 and 160 /jmto accurately estimate the total star formation rate of the Luminous 
Infrared Galaxy population of 9 clusters with mass ~ 10'^ Mq, and 9 groups/poor clusters with mass ~ 5 x 10'^^ Mq. Estimates of the 
field Z(SFR)/M are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the 
universe. 

Results. The field S(SFR)/M increases with redshift up to z ~ 1 and it is constant thereafter. The evolution of the E(SFR)/M-z relation 
in galaxy systems is much faster than in the field. Up to redshift z ~ 0.2, the field has a higher S(SFR)/M than galaxy groups and 
galaxy clusters. At higher redshifts, galaxy groups and the field have similar E(SFR)/M, while massive clusters have significantly 
lower E(SFR)/Mthan both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z ~ 1.6, where the 
group E(SFR)/M lies above the field S(SFR)/M -z relation. We discuss possible interpretations of our results in terms of the processes 
of downsizing, and star-formation quenching. 

Key words. Galaxies: star formation - Galaxies: clusters: general - Galaxies: evolution - Galaxies: starburst 



1. Introduction 

The level of star formation (SF) activity in galaxy systems is 
known to be suppressed relative to the field. Accordingto the 
well known morphology-density relation (Dr essier 19 80 |) and 
the st a r formation rate (SFR)-density rel a tion (|Hashimoto et ahl 
Il998t iLewis et al.1 l2002t iGomez et all l2003h . in the local 
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* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Universe high density regions, like groups and clusters, host 
mostly early type galaxies characterized by a lower level of SF 
activity than field (mostly late-type) galaxies. 

The environmental dependence of galaxy SFR may change 
with redshift, as galaxies in systems undergo significant evolu- 
tion. In higher-redshift clusters, the fraction of blue galaxies is 
higher (the so-called 'Butche r-Oemler' effect lButcher & Oemleii 
119781 11984 IPimbbletll2003l) . and so are the fractions of cluster 
galaxies with spec t ra cha r acterized by young st e llar populations 
(Dr essier & GunnI 119831: iDressler et all 119991: iPoggiand et al] 
119991) and the fraction of galaxies with late-type morphology 
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(iDressleret all 119971: iFasano et alJ 120001: iPostman et al] l2005t 
ISmith et al.l2005h ^Higher-z clusters are also observed to contain 
a higher fraction of infrared (IR) emitting galaxies (the so called 
"IR Butcher-Oemler effect" ISaintonge et al.ll2008L iHaines et aP 
120091: iTemporin et al.ll2009l) . where most of the IR emission is 
powered by SF. 

Given this evolution and the morphology- and SFR-density 
relations observed locally, it is clear that a quenching of the star 
formation activity of galaxies in dense environments is required 
since z ~ 1-2, and that this quenching process must act faster in 
galaxy systems than in the field. When and where this quench- 
ing process takes place, is still a matter of debate. SF quench- 
ing is generally assumed to occur within the cluster environ - 
ment, where processes like ram pressure jGunn & Gottl[T97l . 
cumu lative g alaxy-galaxy hydrodynamic/ gravitational interac - 
tions (iPark & Hwang 200 9). strangulation dLarson et al.lll980l) . 
and galaxy hara ssment dMoore et al.|[T996l) are particularly ef- 
fective. Indeed, iTran et a" ( 20091) observed a higher fraction 
of star forming galaxies in a supergroups at z = 0.37 with 
respect to clusters at the same redshift. It has however also 
been claimed that SF quenching of cluster galaxies occur in 
low-mass groups or large-scale filame nts prior to cluster as- 



semb l y (the so-called "pre-processing" [ Zabludoff & Mulchae 



19981 iKodama et al.ll2o6ll: iFadda et arjl2008t iPorter et alJl2008l 



Balogh et alil2oT l^The accretion process of groups on to clus 



ters c an itself lead to star-formation quenching (Poggi anti et alj 
l2004l) . perhaps caused by rapid gas c onsumption due to a 
sudden enhancement of the SF act ivity dMiller & OwenI [20031 : 
ICoia et al. I I2OO5I: lFerrari"etani200l . It is also not clear whether 
environment plays any role a t all in the quenching process. 
According to lPeng et al.l (|2010|) . mass quenching is the dominant 
quenching process for massive galaxies, which generally reside 
in massive halos like groups and clusters. Since mass quench- 
ing should occur when a galaxy reaches a limiting mass, more 
massive galaxies, which reside in more massive halos, should be 
quenched earlier than less massive ones, having reached earlier 
the limiting mass for quenching. This creates an environment de- 
pendent quenching of galaxy SF, which is driven by some inter- 
nal, rather than external, process, such as AGN feedback. This is 
supported observationally by the analysis of the star-formation 
histories of galaxies in the Virgo cluster region ( iGavazzi et alj 
I2002h . galaxies of higher H-band luminosities being character- 
ized by shorter timescales of SF. Additional support to this sce- 
nario comes fro m the analysis of chemical abundances in ellip- 
tical galaxies bv lPipino & Matteuccil (l2004l) . 

Another way of looking at the evolution of the SF activity 
in galaxy systems is to consider a global quantity such as the 
star formation rate per unit of halo mass, that is the sum of the 
SFRs of all the galaxy in a system, divided by the system to- 
tal m ass, S(SFR)/M. A ccordi ng to recent result s dKodama et alJ 
2004tlFinn et al.ll2004l2005tlGeach et al.ll2006l iBai et al.ll2009t 



Chung et alj 120101: iKovama et alj 120101: faavashi et all 1201 Ih 
the evolution of the cluster S(SFR)/Mis rapid, with a redshift 
dependence 2(SFR)/M oc (1 + z)", with a ^ 5 - 7. The clus- 
ter 2(SFR)/M evolves faster up to z ~ 1 than the fraction of 
IR-emitting galaxies (the IR Butcher-Oemler effect) because IR- 
emitting galaxies are not only more numerous in higher-z clus- 
ters but also more IR lu minous, given the evolution of the cluster 
IR luminosity function (iBai et alJl2009l) . In addition, there are 
some indications that 2(SFR)/M anti-c orrelates with halo mass 
in galaxy systems at similar redshift (Fi nn et all |2004 120051 : 
LHomeier et al. 2005; Koyama et al. 2010). 

The quantity i;(SFR)/M thus appears to be a powerful tool 
for analyzing the rapid evolution of galaxies in systems of dif- 



ferent mass. To determine S(SFR)/M one needs to estimate the 
system dynamical mass and the SFRs of its member galaxies. 
Robust SFR estimates can be obtained from galaxy IR luminosi- 
ties, LiR, and these estimates are considered to be more reliable 
than the SFR based on optical emission-lines such as [Oil] and 
Ha (iGeach et al.1120061) . 

If one considers only the systems where the total SFR has 
been derived from IR observations, the £(SFR) /M-z relation 
is cur rently based on ~ 15 galaxy clusters only ( iBiviano et alj 
l201lh . all with masse£] M2()o ^ 10'"* Mq and z < 1. Masses for 
these systems have been derived in different ways (from X-ray 
or lensing data, or from kinematical analyses of the population 
of cluster galaxies), and at different limiting radii. 

The aim of this paper is to extend the study of the 
S(SFR)/M-z relation to systems of lower masses than clusters, 
and to higher z, with a homogeneous estimate of the system 
masses and total SFRs. Specifically, we extend the mass range 
down to the group mass regime, and to the typical dark matter 
halo mass of field galaxies (lO'^ - 10'^ Mq). We also extend the 
redshift range up to z ~ 1.6. We use deep far-infrared PACS 
data at 100 and 160 //mto accurately estimate the total SFRs of 
18 systems with masses in the range 2 x 10'^* - 3 x IO'^'Mq, 
down to LiR= IO^Lq, i.e. the luminosity that characterizes the 
so-called Luminous IR Galaxies (LIRGs). This limit is set by the 
faintest luminosity observable with our data at z~l.6. We use 
nearly complete galaxy spectroscopic samples to measure the 
system masses from their kinematics. To extend the Z(SFR)/M- 
z relation to the field regime, i.e. to dark matter halo masses 
typical of field galaxi es, we use published global SFR densities 
at different redshifts dMadau et alj|1998l: lOruppioni et alj|201 U 
Magnelli et al. 20ll|). 

In Sect. 12.11 and Sect. 12.21 we describe our sample and our 
spectroscopic and photometric data-set, in Sect. 12.31 we de- 
scribe how we define system membership of the galaxies with 
available redshifts, and derive the system dynamical masses, in 
Sect. l2.4l we determine the total SFRs of the systems. The result- 
ing S(SFR)/M-z relations are presented in Sect. [3] for clusters, 
groups (Sect, [ni l, and the field (Sect. 1X2] ). We summarize and 
discuss our results in Sect.|4] 

We adopt Ho = 70 km s"' Mpc"', 0„, = 0.3, Q.a = 0.7 
throughout this paper We adopt a Salpeter IMF when calculating 
stellar masses and SFRs. 



2. The data-set 

2. 1 . The sample of galaxy systems 

Our sample comprises 9 X-ray detected systems at 0.15 < 
z < 0.85, all obser ved in the PACS Evolutionary Probe PACS 
GT Program (PEP. iLutz et alJl201 ih . The systems are listed in 
Table [1] In addition, we use the X-ray detected group and clus- 
ter catalog of the COSMOS field (Finoguenov et al., in prep.) 
to identify 27 systems at z < 0.8 with at least 10 spectro- 
scopic members each, the minimum requir ed for an acceptabl e 
mass estimate from the system kinematics (iGirardi et alj|1993h . 
The COSMOS field is observed as part of the PEP program. 
We further include in our sample 4 systems in the GOODS 
North and South fields at 0.7 < z < 1.6 observed as part 
of the PACS GOODS-Herschel Program (Elbaz et al. 2011). 



' The cluster virial mass, M200, is the mass contained within a sphere 
of radius r2oo> which is the radius within which the enclosed average 
mass density is 200 times the critical density. The two quantities are 
related by the expression M200 = 100H(z)'r2p(,/G, where G is the grav- 
itational constant and H(z) the Hubble constant at the system redshift. 
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Of these sy stems, two are X-ray detected, one at z - 0.73 in 
GOO DS-S (iGilli et al.l2003HLe Fevre et al.l2004tlPopesso et al 



20091) . another one at z = 1.02 in GOODS-N (lElbaz et al 



2007h . one system lies close to the Chandra CCD chip gap 
and it is not X-r ay detected (GOODS-N system at z = 0.85 
iBauer et"al ]l200l .' and the highest redshift system (at z = 1 .612) 
is on ly m arginally X-ray detected in the Chandra GOODS-S 
jKurk et al. 



cti 

>. 



10' 



map 



.2008.). 



Finally, we include in our sample the 

Bullet cluster (iBarrena et al.l 120021: iMarkevitch et all |2004 . a 



well-known massive cluster undergoing a violent collision with 
an infalling group, which we use as a test-case for the depen- 
dence of the 2(SFR)/Mon the cluster dynamical state. PACS 
data for the Bullet cluster are not public yet. For this cluster 
we use the Lir base d on Spitzer MIPS 24 fiva data published by 
IChung et al.l (l2010h . The considered LlRGs members lie within 
1.7 Mpc from the cluster center This is close to the r20() ra- 
dius, ~ 2 Mpc, estim ated from the kinematics of cluster galaxies 
(IBarrena et al.ll2002l) . 

We split our sample in groups and clusters depending on 
their mass estimates (see next section), groups/poor clusters with 
10'^^ Mo < M200 < 3 X IO'^^Mq and clusters with M200 > 
3 X 10'"^ Mq. The mass value separating groups/poor clusters 
from clusters ap proximately corres ponds to the richness class 
in the catalog of lAbell et al.l (Il989l) (using the relations between 
mass and velocity dispersion and between velocity dispersion 
and richness class for galaxy clusters; see iGu'ardi et al..,1993t 
iBiviano et al I I2OO6I) . 



2.2. Spectral and photometric data 



10' 



10' 




M200 (from CT„ 



Fig. 1. Comparison of system mass estimates obtained from X- 
ray data with those obtained from system velocity dispersions. 



men ted with MI PS 24 ^m data, \y hen available, and PACS 100 
and 160 //mdata ( iLutz et al.ll2011h . 

The data reduction, source extraction, and a detailed expla- 
nation o f the depths of the PAC S maps used in this paper can be 
found in lLutz etalj ( 1201 lb and lElbaz et all (1201 ll) . for the PEP 
and GOODS -Herschel surveys, respectively. 



We take the data for the COSMOS field from the multi- 2.3. Membership and mass estimates of galaxy systems 



wavelength and IRAC-selected catalog of lllbert et al.l (120 IC ) 
complemented with public Spitzer 24 yum dLe Floc'h et al.ll200S ; 
Sanders et al. 2007) and PACS 100 and 160 ^mdata dLutz et all 
201 ll) . The association between 24 yumand PACS sources 



with their optical counterparts is done via a maximu m likeli- 
hood method (see iLutz et alj 120111: iBerta et alj 1201 Ol for de- 
tails). The photometric sources were cross-matched in coor- 
dinates with the sources for which a high-confidence spec- 
troscopic redshift is available. For this purpose we use the 
public catalogs of spectroscopic redshifts complemented with 
other unpublished data. This catalog includes redshifts from 
either SDSS or the public zCO SMOS -bright data acquired 
using VLT/VIMOS dLillv et al.l l2007l 12009) complemented 
with Keck/DEIMOS (Pis: Sco yille, Capak Salva to, Sanders, 
Kartaltepe), Magella n/IMACS (iTrump et all l2007h . and MMT 
(iPrescott et al.ll2006l) spectroscopic redshifts. 

Similar multi-wavelength (from UV to 24 /im) data for 
the GOODS-S field is provided by the MUSIC catalog 
dGrazian et al.l 120061 : ISantini et alJ l2009l) and complemented 
with PACS data at 70, 100 and 160 jumand spectroscopic red- 
shifts from the PEP catalog (se e Lutz et al. 201 1 fo r more details) 
and from the GMASS survey dCimatti et al.ll2008l) . 

We take the data for the GOODS-N field fro m the multi- 
wave length catalog prepared by the PEP Team dBerta et alJ 
l20Toh including UV to MIPS 24 data, PACS 100 and 160 
^ mdata, complement ed with the spectroscopic redshift catalog 
of lBarger et aD d2008l) . 

We collect publicly available spectroscopic data for all clus- 
ters from the NED. For CL 0024-H1 7 and MS 0451-03 we als o 
collect multi-wavelength data from iMoran et all d2005L l2007l) . 
Multi-wavelength photometric data are not available for the 
other clusters. The references for each cluster spectroscopic 
sample are listed in Table [Tj Each cluster sample is comple- 



To define the galaxies that are members of the systems describe d 
in Sect. 12.11 we adopt the algorithm of iMamon et alJ d2010l) . 
which is based on the modeling of the mass and anisotropy pro- 
files of cluster-sized halos extracted from a cosmological nu- 
merical simulation. This algorithm is more effective than tradi- 
tional approaches (e.g. iYahil & Vidal 1977 ) in rejecting interlop- 
ers, while still preserving cluster members. The system member- 
ship selection depends on the location of galaxies in the system- 
centric distance - rest-frame velocitjQ diagram. The peaks of the 
X-ray surface brightness are adopted as centers of the X-ray de- 
tected systems. The galaxy number density maxima (estimated 
using an adaptive kernel technique) are adopted as centers for 
the other systems. The interloper rejection procedure is iterated 
until convergence. 

As an outcome of the procedure of iMamon et al. I (l20Toh . 
the syst em dynamical ma sses are obtained from biweight es- 
timates (iBeers et al.lll990l) of t he system velocity dispers ions 
along the line-of-sight, cTy, as in iMauduit & MamonI (l2007l see 
Appendix A). We use the cTy-based mass estimates, rather than 
the more traditional virial mass estimates because the o-y-based 
mass estimates are robust against problems of spatial incom- 
plet eness, unlike the mo re traditional virial theorem estimator 
(see iBiviano et al" ] l2006l) . so we do not have to worry about the 
possibility that in some systems of our sample the spectroscopic 
completeness is different in different regions. 

The uncertainties in the system mass estimates a re estimated 
from the uncertainties in the biweight estimate of cry dBeers et al.l 
119901) via the propagation of error analysis. 



- The galaxy rest-frame velocities with respect to the system 
mean velocity are obtaine d by the usual relation v = c(z - z)/(l + z) 
dHarrison & Noonaiilll979l) . where z is the system me an redshift, deter- 
mined with the biweight estimator dBeers et aDll990h . 
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name 


z 


(Mpc) 


M200 
(lO'^Mo) 




NuRG 




ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Cluster sample 


A2218 


0.175 


2.73 


27.6 


107 


1 


75% 


1 


A2219 


0.226 


2.375 


19 


126 


2 


68% 


2 


A2390 


0.228 


2.533 


23 


259 


3 


69% 


3 


MS 1358+62 


0.329 


1.711 


8.2 


286 


1 


73% 


4 


COSMOS.CLl" 


0.351 


1.2 


4.8 


36 


6 


48% 


5 


MS 0451-03 


0.539 


1.966 


15.5 


342 


8 


80% 


6 


COSMOS.CL2" 


0.78 


1.18 


3.88 


36 


6 


48% 


5 


MS 1054-03 


0.830 


1.521 


10.1 


134 


5 


100% 


7 


RX J0152-13 


0.838 


1.657 


13.3 


228 


11 


96% 


8 


Group/poor cluster sample 


COSMOS.GRl" 


0.1 


0.61 


0.26 


47 


1 


56% 


5 


COSMOS.GR2'' 


0.22 


0.62 


0.31 


77 


1 


66% 


5 


COSMOS.GR3'' 


0.35 


0.68 


0.31 


81 


4 


48% 


5 


CL 0024+17 


0.395 


0.96 


1.5 


396 


10 


92% 


6,9 


COSMOS.GR4'' 


0.70 


0.69 


0.31 


34 


10 


47% 


5 


GOODS - S.Ol 


0.735 


0.51 


0.85 


188 


3 


88% 


10 


GOODS - N.Ol 


0.85 


0.51 


0.9 


112 


3 


90% 


11 


GOODS - N.02 


1.02 


0.73 


1.03 


130 


4 


79% 


11 


GOODS - S.02 


1.61 


0.52 


0.9 


76 


6 


68% 


10 


Merging system 


Bullet cluster 


0.297 


1.7 


9.5 


43 


5 
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Table 1. The table lists the following main properties of the cluster and group samples: column (1): name; column (2): redshift; 
column (3): radius raoo in Mpc; column (4): mass M200 in 10'"* M©; column (5): number of members used in the estimate of 
M200; column (6): number of LIRG members w ithin T200; column (7): sp ectroscopic co mpleteness of the LIRG s ample within r2oo ; 
c olumn (8): referenc es for spectroscopic data: 1 iLe Borgne etaP d 19921); 2 SPSS DR7 (|Abazaiian et al."2009); 3'Yee et al.' ('1996|); 
4lFisheretal.[(ll998h ; 5 COSMOS spectroscopic sample, see text; 6 lMoran et al.l(l2005l);^ Tran et al. (2007); 8 Pat el et al. (2009); 9 
iMoran et ain2007[); 10 P EP- GOODS-S spectrop hotometric catalog iPopesso et all (120111) ; 1 1 PEP-GOODS -N spectrophotometric 



catalog lBerta et al 



(120101) : 12 lChunget all (120101) . 



a Composite of COSMOS clusters. The quantities reported for this system are the mean of the contributing systems, 
b Composite of COSMOS groups. The quantities reported for this system are the mean of the contributing systems. 



As an alternative mass estimate we consider that obtained 
by assuming hydrostatic equilibrium of the X-ray-emitting intra- 
cluster or intra-group gas. The two mass-estimates are compared 
in Fig.[T]for those systems which have both estimates available. 
The scatter is large, but there is no significant bias, on average, 
between the two mass estimates. There is only one system for 
which the two mass estimates are in disagreement by more than 
3cr. The discrepancy is probably due the low number of spec- 
troscopic members (9) available for this system. We exclude it 
from our sample for further analyses. 



2.4. Star formation rates 

We fit the spectral energy distributions (SEDs) of all sys- 
tem members in our sample which have multiwavelength data 
available from the UV to the PACS wav elengths, using a set 
of empiri cal templates of local objects dPoUetta et al. l2007t 
iGruppioiii et al. 2010) that reproduce most of the observed 
galaxy SEDs from th e optical to the FIR wavelengths (see 
iRodighiero et al.l2010l) . We then determine the galaxy Lir by in- 
tegrating the best- fit SED models from 8 to 1000 /^m. 

We estimate the Lir of those system members for which mul- 
tiwavelength data are not available, by fitting their PACS data 
with SED models from Polletta et al. (2007), to extract the vLy 
at 100 yumrest frame, and by using the Ljr -vLy relation, as ex- 
plained in Popesso et al. (in prep.). 



Some of our member galaxies do not have PACS data ei- 
ther due to the incompleteness of the PACS photometric cat- 
alog between the 3 and Scr levels, or because the PACS cat- 
alog are not deep enough to reach the required 1O"L0 lim- 
its at the redshift of the system. We determine the Lir of these 
galaxies from their MIPS 24 ^m flux densities by using the scal- 
ing relation of lCharv&Eiba3(l200Ih . corrected as ex plained in 
Popesso et al. (in prep.) to avoid overestimations (iElbaz et al.l 
I2OI UlNordon et al.l201 Ih . We also re-determine in the same way 
the Lir of t he member galaxies of the Bullet cluster from the 24 
//mdata of lChung etall (l2010h . In total, only 8 galaxies (9% of 
our member LIRGs) do not have their Lir estimated from PACS 
data: 2 in MS 0451-03, 4 in RX J0152-13 and 2 in the z ~ 1.6 
GOODS-S group. 

The main advantage of using PACS data in the determina- 
tion of the Lir is to remove possible contamination by the emis- 
sion of active galactic nuclei (AGN), since most of the rest- 
frame far-IR emission comes from the host galaxy (iNetzer et al.l 
I2OO7I: iLutz et al.ll20Toh . Thus, we do not need to exclude AGN 
hosts from ou r sample as was instead done in previous similar 
analy sis (e.g. iGeach et al.ll2006t iBai etal.ll2009t iBiviano et all 
1201 ih . AGN contamination is a concern, however, when esti- 
mating the SFR from 24 fim data. We do exclude the AGN 
hosts from the evaluation of the total SFR of the Bullet clus- 
ter, using the AGN classification of Chung et a l. (2010), since 
PACS data are not yet publicly available for this cluster. We 
also cross-correlated the system members, whose Lir is based 
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on 24 ^imdata only , with publicly avail a ble AGN catalog s 
(iJohnson et alJ |2003[ iDemarco et all 120051: iBrusa et al.l l2009h . 
None of those members turn out to be an AGN. 

Finally, we estimate the galaxy SFRs from their Lir via the 
law of Kennicutt (1999). The LIRG LiR-limit corresponds to a 
Umiting SFR of 17 Mq yr^'. 

We estimate Z(SFR), the total SFR of each system, by sum- 
ming up the SFRs of the LIRGs that are members of the system. 
We then correct 2(SFR) for spectroscopic incompletetness in the 
following way. For all systems where the PACS data reach the 
1O"L0 limit at the system redshift, we take as a reference the 
PACS 100 /im catalog. We plot the 100 yum flux density vs. Lir of 
system member galaxies in narrow redshift bins to estimate the 
100 /imflux density that corresponds to the LIRG LjRat each 
redshift, fuRcCz). The spectroscopic incompleteness is defined 
as the fraction of PACS sources with z, among all those with flux 
> fLiRG(z) within r20(). The inverse of this fraction is the correc- 
tion factor that we apply to the S(SFR) estimates to correct for 
incompleteness. For the systems which contain member galax- 
ies for which the LiR-estimates are based on MIPS 24 jumflux 
densities, we use the same approach by using the 24 jim catalog 
as reference. 

To test the reliability of this method, we consider a 
system with 100% spectrosco pic completenes s down to the 
LIR Lir level, MS 1054-03 dBai et al.l l2009h . We perform 
500000 random extractions from the luminosity function of 
MS 1054-03 of a number of LIRGs corresponding to a given 
spectroscopic completeness. We repeat this test for all the spec- 
troscopic coverage values listed in Table 1 . We then estimate the 
uncertainty due to the incompleteness correction from the dis- 
persion of the difference between the recovered E(SFR) and the 
real value. The uncertainty due to this correction is ^ 20% in the 
range of spectroscopic completeness of our sample. 

The total uncertainty of the 2(SFR) estimates is determined 
from the propagation of error anal ysis, by consider ing a 10% 
uncertainty in the Ljr estimates (see lLutz et al.ll201ll for further 
details), and the 20% uncertainty due to the completeness cor- 
rection. 

The spectroscopic coverage of the COSMOS field is much 
lower than those of the clusters and GOODS fields and spatially 
not uniform. Thus, to reduce the error bars of the 2(SFR), we 
combine the COSMOS groups and, separately, the COSMOS 
clusters in four, and, respectively, two redshift bins, see Table [T] 
In each z-bin we then define the mean E(SFR), M2()o, and spec- 
troscopic incompleteness, separately for the groups and the clus- 
ters contributing to the sample in that bin. 

Counting the composite COSMOS groups and clusters as in- 
dividual systems, our final sample comprises 9 massive clusters, 
9 groups/poor clusters, and the Bullet cluster (see Table[T]). 



3. Results: the i:(SFR)/M-z relation 

3. 1 . Clusters versus groups 

The S(SFR)/M-redshift relation is shown in Fig. |2] Black, 
and magenta symbols show the Z(SFR)/M- redshift relation 
for the clusters and, respectively, the groups. For both sam- 
ples there is evidence for a significant 2(SFR)/Mvs. z corre- 
lation, 99% significant according to a Spearman rank correla- 
tion test. On the other hand, we do not find evidence for sig- 
nificant S(SFR)/M-M2()o or M2()o-z correlations within the group 
and cluster samples separately. Hence the observed Z(SFR)/M- 
z relation within each sample must be interpreted as a genuine 
redshift evolution of the E(SFR)/M of galaxy systems. 
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Fig. 2. E(SFR)/M- redshift relation for clusters (black symbols) 
and groups (magenta symbols). Square symbols, triangles, and 
dots identify, respectively, the COSMOS composite systems, the 
GOODS systems, and the remaining systems. The red star iden- 
tifies the Bullet cluster. The black solid line shows the best fit 
S(SFR)/M-z relation for the cluster sample, excluding the Bullet 
cluster. The relation fitted to the data is of the type S(SFR)/M oc 
f. The dashed l ine shows the relation, S(SFR)/M oc (1 + zY , of 
iBai et al.l (|2009|) rescaled to the LIRG LiR-regime and to the r2oo 
region. The magenta solid lines shows the best fit for the sample 
of groups and poor clusters. 



Groups appear to be characterized by higher 
2(SFR)/M values than clusters, at all redshifts, i.e. they 
show a higher SF activity than massive clusters. In this sense, 
when the sample of groups and clusters are considered together, 
S(SFR)/M does anticorrel ate with the sys t em M200, as already 
found in previous studies dFinn et al ]|2004 120051: 1 iHomeier et al.l 



l2005llKovama et al.ll20Tol) . 

The solid curves in Fig.|2]represent best-fit models to the ob- 
served S(SFR)/M-z relations, I(SFR)/M = (66 + 23) x zi "±0-36 
for the cluster sample and S(SFR)/M = (213+44) x z'-33±0-34 for 
the group sample. In Fig. |2]we also sho w as a dashed lin e the 
2;(SFR)/Moc (1 + z)5^ model proposed bv lBai et al.l ( |2009|) . We 
rescale their relation, which is built for S(SFR) computed down 
to a SFR limit of 2 M© yr ', to ou r adopted LIRG limit, by us- 
ing the IR luminosity function of iBai et"ar l (l2009l) . Moreover, 
the 2(SFR) in Bai et al.'s relation is evaluated within a radius 
of rsoo, the radius corresponding to an overdensity equal to 500 
times the c ritical density. To co nvert to our adopted radius r2oo, 
we adopt a iNavarro et al.l (Il997l) model, wich concentration c=5 
to scale for the number of galaxies inside the two radii and de- 
rive the r5oo/r2oo ratio, and then use Fig. 7 in iBai et all (|2009|) 
to account for the different fraction of IR-emitting galaxies at 
rgoo and r2oo- The modified relation of _Bai et al. (20()9j) appears 
to fail to predict th e rapid evolution of c luster S(SFR)/M, as al- 
ready suggested bv lBiviano et al ] (1201 11) . 
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3.2. Structures versus field 

In Fig. |3] we compare the Z(SFR)/M-z relation of galaxy sys- 
tems with the corresponding relation for field galaxies (light 
blue shaded region and blue dashed line). The 2(SFR)/M-z re- 
lation of field galaxies is obtained b y dividing the observe d Star 
Formation Rate Density (SFRD) of|Magneni etaD (1201 iL trian- 
gles) and the modeled SFRD of iGruppioni et al.l (1201 iL dashed 
line), by the mean comoving density of the universe (Q,„ x pc 
where Q„, = 0.3 and pc is the critical density of the Universe). 
Both SFRD have been evaluated only down to the SFR corre- 
sponding to the LIRG Lir , via the Kennicutt relation. 

The field SFRD has been estimated in large comoving vol- 
umes that include galaxy systems, voids, and isolated galaxies, 
and is thus representative of the general field g alaxy popula- 
tion. A ccording to the dark halo mass function of Ijenkins et al.l 
(l200lh . halos of 10" - 10'- Mq give the main contribution to the 
dark matter budget at all redshifts. Thus, the Z(SFR)/M of Fig. [3] 
can be considered as an effective estimate of the S(SFR)/M of 
galaxy-sized dark matter halos. 

Observed and modeled Z(SFR)/M are in very good agree- 
ment within the error bars. They increase from z = to z ~ 1 
where they reach a plateau. 

The field Z(SFR)/M-z relation lies above both the group and 
the cluster relations at z ^ 0.2. The field E(SFR)/Mis more 
than an order of magnitude higher than those of Abell 22 1 8 at 
z = 0.175 and of the composite COSMOS group at z = 0.1. 
At higher redshifts, group and field galaxy halos show compara- 
ble 2;(SFR)/M. Possibly, the group Z(SFR)/MHe above the field 
relation at z ~ 1.6, consiste ntly with a reversa l of the SFR- 
density rela tion observed by Elbazetal.1 (l2007h . ICooper et alJ 
(120081) . and iPopesso et al] (1201 ij). However, this conclusion is 
based on a single group S(SFR)/M determination. Clusters have 
lower S(SFR)/M than field galaxy halos at all redshifts up to the 
last measured point at z = 0.85. A blind extrapolation of the 
best-fit cluster E(SFR)/M-z relation would suggest that clusters 
should display a higher SFR per unit mass than field galaxies at 
at z ^ 2. However, such an extrapolation is extremely uncertain 
with the present data and massive systems as the ones considered 
here should not even exist at such redshift according to ACDM 
hierarchical model predictions. 



4. Discussion and conclusions 

We find significant evidence for an evolution of the 
S(SFR)/M with z in clusters, groups, and the field. Both the clus- 
ter and the group 2(SFR)/M increase monotonically with z and 
the field 2(SFR)/M reaches a maximum at z ^ 1 . 

The cluster S(SFR)/M evolution we find, appears to be faster, 
at least in the range < z < 1, compar e d to the evolutionary re - 
lation suggested bvlGeach et al.1 (l2006l) : lBai et alj (f2007ll2009h : 
lKovamaetai](l2010ir This might be due to our higher-LiR cut, 
if LIRGs evolve faster than less luminous IR-emi tting galaxies. 
A fast er evolution has already been suggested by iBiviano et alj 
(1201 Ih . but only in the z 5 0.5 range, while in our data there is 
no sign of a flattening of the S(SFR)/M-z relation for clusters 
above that redshift. 

We confirm earlier ind ications dFinn et alJ l2005t iBai et alJ 
l2007t iKovama et al.ll20Tol) that S(SFR)/Mis lower for systems 
of higher mass, and we show for the first time that this is true 
at all redshifts from z ~ to z ~ 0.9. This is true when we 
compare clusters to groups and the field. On the other hand, the 
group S(SFR)/Mlies below the field relation only at z $ 0.2. 
At higher redshifts, the group and the field 2(SFR)/M become 
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Fig. 3. The field 2(SFR)/M-redshift rela tion fromlMagnelli et al.l 
(1201 ih (light blue shaded region) and IGruppioni et alJ ( 201 1 ) 
(dashed blue line). The shading and error bars represent Icr con- 
fidence levels Other symbols/lines have the same meaning of the 
symbols as in Fig.|2] 



comparable. There is a hint of a reversal of the SFR activity vs. 
environment at z ~ 1 .6, where the group E(SFR)/M lies above 
the field S(SFR)/M-z relation. 

Our results appear to support a scenario in which the quench- 
ing of SF occurs earlier in galaxies embedded in more mas- 
sive halos, i.e. first in clusters, then in groups, and finally in 
the field. This would be consistent with a "halo downsizing" 
effect, wher eby massive halos e volve more rapidly than low 
mass halos dNeistein et al.ll2006l) . Although in a narrow range 
of halo masses, the "halo downsizing" effect has already been 
observed in th e stella r-to-total mass ratio vs. redshift relation by 
iFoucaud et all (120101) . 

What is causing the difference between the cluster and 
the group/field E(SFR)/M-z relation? The similarity of the 
field and the groups Z(SFR)/M-z relations suggests that a 
SF quenching effect is taking place mostly after galaxies en- 
ter the cluster environment, and not in groups before they 
merge into more massive structures, as wo uld be predicted 
by the pre-processing scenario ( Zabludoff & MulchaevI 119981 : 
iKodama et al1l200ll:lBalogh et alfcol II) . We can therefore con- 
sider the field and group S(SFR)/M-z relation as indicating the 
normal galaxy evolution, and the cluster 2(SFR)/M -z relation as 
the accelerated evolution experienced by galaxies following their 
accretion to more massive structures, as individuals or in groups. 
However, we point out that other differences may intervene be- 
tween the star forming population of the groups and the field. 
Indeed, in agrement with our results. lTran et al.l (l2009l) find con- 
sistent levels of star formation in the supergroup at z = 0.37 and 
the field population, but they also reveal a much higher fraction 
of early type galaxies among the SF galaxies of the group en- 
vironment with respect to the field. Similar results are obtained 
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bv lHwang&ParkI (120091) . This could indicate a morphological 
transformation happening prior to the SF quenching effect. 

Part of the observed evolution may come from internal, 
rather than external processes. In fact, most of the LIRGs in 
our sample are likely to be rather massive galaxies (as inferred 
from the almost linear, z- dependent, relatio n between galaxy 
SFRs and stellar masses o f lElbaz et al.ll2007h . and according to 
iPeng et all (l2010h massive galaxies evolve mostly because of an 
internally driven process, called 'mass quenching', caused per- 
haps by feedback from active galactic nuclei. But since this pro- 
cess is unlikely to be more efficient in quenching SF of massive 
galaxies in clusters than in other environments, we need to con- 
sider an additional quenching process that would operate only in 
the cluster environment. 

Two obvious candidates for producing this addit ional 
quenching, are ram-pressure s tripping (Gunn & Goti l 19721) and 
starvation ( Larson et al. 19801) . Ram-pressure stripping is a fast 
process dAbadi et al. 19991) and is effective only in galaxy sys- 
tems where the gas density is high, i.e. in clusters. However, 
ram-pressure stripping seems ruled out by the presence of LIRGs 
contributed by the infalling group (the "bullet") in the Bullet 
cluster The group LIRG contribution is probably the reason 
of the position of the Bullet cluster in the 2:(SF R)/M-z plot 
interm ediate between that of clusters and groups dChung et al.l 
I2OIOI) . If rapid ram-pressure stripping was a major star formation 
quenching process, the group LIRGs could have had their star 
formation quenched already, g iven that the cluster-gro up colli- 
sion happened ~ 250 Myr ago (ISpringel & Farrai]l2007h . 

Starvation, caused by the removal of the hot gas halo reser- 
voirs of galaxies, is a more likely candidate given that we see an 
accelerate d evolution (relative to the field) both i n clusters and 
in groups dKawata & Mulchaevl2008T:lBekkil2009l) . The removal 
of galaxy hot gas reservoirs inducing starvation can be caused 
by tidal galaxy-galaxy encounters or by the interaction with 
the intra-cluster/intra-group medium. Starvation should proceed 
more effectively in higher (galaxy or gas) density regions, hence 
it should quench SF earlier in cluster than in group galaxies, as 
we observe. 

In order to better constrain these evolutionary scenarios it 
is important to extend the analysis of the Z(SFR)/M -z relation 
to even higher redshifts, where we might expect to see a rever- 
sal of the environment dependence as we approach the epoch 
when galaxies in dense environments underwent their first ma- 
jor episodes of SF. In the near future we plan to use new Herschel 
PACS observations of 6 poor and rich clusters at 1 .4 < z < 1 .9 
to explore this issue. 
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